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I review what is required to compute transport coefficients in ultra-relativistic, 
weakly coupled gauge theories, at leading order in g, using kinetic theory. Then I 
discuss how the calculation would look in alternative approaches: the 2PI method, 
and direct diagrammatic analysis. I argue that the 2PI method may be a good 
way to derive the kinetic theory, but is not very useful directly (in a gauge theory). 
The diagrammatic approach is almost hopeless. 


1. Introduction 

A “big question” in particle physics at extreme energy densities (Strong 
and Electroweak Matter, the topic of this conference), is how nonequilib¬ 
rium systems behave. Here I will only talk about weak coupling, plus some 
conditions on the departure from equilibrium and the measurables of in¬ 
terest, which allow application of techniques grounded in perturbative field 
theory. (Some very interesting questions are excluded by this condition; 
but the problems included are interesting and challenging.) In particular, 
I am going to talk about transport coefficients, even though they represent 
the “closest to equilibrium” nonequilibrium systems we could consider. 

Transport coefficients tell about equilibration in systems which are ho¬ 
mogeneous on fairly large scales, and are locally fairly close to equilibrium. 
The traditional approach to dealing with such a system at weak coupling is 
kinetic theory, that is, Boltzmann equations. We have recently presented 1 a 
Boltzmann equation which is adequate for treating transport coefficients in 
a nonabelian gauge theory at leading order in the coupling g (see also Larry 
Yaffe’s talk in this proceeding). Its applicability extends well beyond the 
near-equilibrium situation needed for transport, though the homogeneity 
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requirements for the equation to apply are stronger than you might naively 
guess. 

But there has been much interest in addressing transport using other 
tools, particularly the 2PI formalism and direct diagrammatic analysis. I 
will argue that, if the goal is to perform a calculation in a gauge theory, 
which is complete at leading order in the gauge coupling g, then both meth¬ 
ods become very complicated. The 2PI formalism may be a very efficient 
way of deriving the correct Boltzmann equation, but direct application 
without gradient expansion-for instance, numerically, as has recently been 
done in scalar field models 2,3,4 -will encounter substantial difficulties. The 
diagrammatic approach at leading order appears to be so complicated that 
I see little reason to pursue this method further. 


2. Boltzmann approach to transport coefficients 

When a plasma is sufficiently homogeneous, the large scale flow is accurately 
described by ideal hydrodynamics. Under ideal hydrodynamic behavior, 
entropy is conserved; particle number is often approximately conserved. In 
the heavy ion context, this means that almost all the available energy of 
the system goes into bulk flow. 

Transport coefficients give the first corrections to this behavior, when 
the system is large but not infinitely large. They cause entropy generation, 
and tell how close to local thermal equilibrium the system remains. In the 
context of heavy ion collisions, knowing the transport coefficients could tell 
us when the hydrodynamic approximation has broken down, which would 
help hydro simulators to know when to impose chemical and kinetic freeze- 
out. In early universe physics, diffusion coefficients are useful in baryogene- 
sis calculations, and viscosity and diffusion coefficients can indicate how far 
around electroweak bubble walls the departure from equilibrium extends. 
For the “pure” theory of thermal field theory, transport coefficients are a 
useful object of study because their definition is very theoretically clean. If 
we want to develop the tools to study nonequilibrium problems in a con¬ 
trolled way, those tools should be able to tell us the values of transport 
coefficients. So being able to compute them is the most basic step towards 
more general nonequilibrium calculations. 

In the Boltzmann approach, one shows (or assumes, or hopes) that the 
plasma is well described as a bath of long lived quasiparticles, which un¬ 
dergo occasional and relatively localized scatterings. That is, one requires 
that the interesting observables are dominated by the two point function, 
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that the spectral weight for the two point function is tightly peaked, that 
the evolution of the two point function can be accurately written in terms 
of scatterings with the other quasiparticles, and that these scattering pro¬ 
cesses are local on the scale of the spatial variation of the system. 

Given these conditions, one can write down a Boltzmann equation, 

d t f(p , *>*) + “?• d *f(P> t) + F 9 p f(p, x , t) = —C[f] , (1) 

with C the collision operator. Applying this to a situation where / is 
spatially inhomogeneous, or where the force term represents an electrical 
field, allows calculation of transport coefficients. In particular, when / 
exhibits shear flow, then computing the traceless part of the stress tensor 
gives the shear viscosity, 

Tij[f] - i 5ijT kk [f} = r] (diVj + djVi - ^SijdkVkj . (2) 

Similarly, the diffusion constant is determined by the size of a conserved 
current, when the number density varies in space. In each case, one must 
solve for / (at linearized order in the departure from local equilibrium) in 
the Boltzmann equation. This is hard because it is an integral equation; 
the operator C[f] contains / at momenta other than p. It is especially hard 
because, at leading order, C turns out to be quite nasty! 

The collision operator C should include “all collision processes which 
are important at leading order.” At leading log order, there are only a 
few processes; t channel boson exchange, Compton scattering, and an an¬ 
nihilation process related to Compton by crossing. Further, simplifying 
approximations can be made to each collision term, which in fact makes 
the collision operator a differential operator, not an integral one 1 . But at 
leading order, in fact at next-to-leading-logarithm, one needs not only all 
2 2 processes, but certain inelastic LPM suppressed splitting processes 5 . 

For the full details see 5,6,7 . To show how bad the situation is, I will write C 
in some detail, though not enough to allow you to evaluate it: The 2 <-> 2 
piece is 

Cr 2 [/] = t4-E / \Mt d {p,k-,p',k')\\2n)H^{P+K-P'-K') 

° bed Jk P' k ' 

x{fa( P ) fb(k ) [1 ±fc( P ')} [1 ±fd(k')} - [(p, k ) - (p\ k ')]} . (3) 

The 2 <-> 2 matrix elements should be summed over incoming and outgoing 
spins and colors; literature values, summed only over outgoing states, are 
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available 8 . But in addition, HTL corrections are needed on some inter¬ 
nal lines, which means that the form of certain matrix elements is more 
complicated than in vacuum 5 . But the real problem is the inelastic LPM 
suppressed splitting processes: 

C‘a~ 2 if] = I dp ' dk ' ‘HIpI - p' ~ k "> 7b c (p;p' P, k' p) 

Z \P\ V a b c JO 

x{/a(p)[l ± fb{p' P)][ 1 ± fc(k’p)\ -[/<“► (1±/)]} 

+ hit /L f dkdp ' ^(IpI + k ~p') i c ab(p'PiPi kp) 
\p\ u *xz J ° 

x{fa(p)fb{kp)[ 1 ± fc{p' P)] ~ [/ ^ (1±/)]} • (4) 


The awfulness is hiding in J bc (-..); evaluating it requires solving an integral 
equation. For instance, for q <-> qg processes, it is 

7 q q ipn-,p'n, kh) = P + P 7q(p,p', k ), 


^(P,p',k) = 


p' 2 p 2 k 3 q 
dpCpas f d 2 h 


2(2 it) 3 J (27r) 

2h = i5EF™{h)+g 2 T 


2 2ft ReF ™{h;p,p',k), 


f d2( l -L1 

rj__ 

_ 1 ^ 

/ (27r) 2 1 

Ui 

q± + to d/ 


x \{C a -\Cx) F?(h)-Fi(h-kq ± ) 




F”(h) - F*(h+pq ± ) 


+iC A F?(h)-F?(h-p?q ± ) 


SE = 


m. 


eff.g 


to ; 


eff ,q 


m 


eff ,q 


+ 


h z 


( 5 ) 


2 k 2 p' 2p 2pkp' 

That is, there is an integral equation just to determine the coefficient in 
one of the collision terms, which itself appears in an integral equation for 
the expression that we want. 


3. Alternate approaches 

For some reason, people don’t like the Boltzmann approach. The obvious 
reason is that it involves solution of an integral equation (in fact, at leading 
order, an integral equation within an integral equation). This objection 
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is probably not well grounded, because other approaches presumably re¬ 
duce to the same integral equation, in one guise or another. But there are 
some deeper, more philosophical reasons people don’t like the Boltzmann 
approach. 

One problem is that people don’t trust it. Where is the careful deriva¬ 
tion? How do you know you have not missed some necessary term in the 
collision integral? What about quantum and high gradient effects? What 
do you do when the particle widths are large? Another problem people have 
is that the Boltzmann approach is old fashioned. Boltzmann was playing 
with such equations a hundred years ago. Surely there is a more powerful 
modern approach. 


3 . 1 . The 2PI approach 

One approach much advocated in the recent literature 9,10 ’ 2,3 ’ 4 (see also 
contributions to this proceedings from Berges and from Mottola) is the 
2PI approach. The idea, in summary, is to solve for the evolution of the 
2-point function by variationally minimizing a functional 12(A) of the 2- 
point function A. This functional looks like a logarithm of the 2-point 
function, plus a term involving the propagator and self-energy, plus <f>, the 
sum of 2 particle irreducible (2PI) bubble diagrams; inside the functional 
the 2-point function is treated as arbitrary. If the complete set of 2PI 
diagrams is included, the method is exact; but generally, to get anywhere it 
is necessary to truncate in some way the sum over 2PI diagrams. There are 
two approaches to using the 2PI formalism which particularly recommend 
themselves: 

(1) Make additional approximations, in particular, make a gradient ex¬ 
pansion and assume slow time development. Similar approximations 
to those Calzetta and Hu used in scalar theory 9 may be the most 
efficient way to derive the Boltzmann equation described in the last 
section. 

(2) Take no prisoners. Directly solve the evolution of the 2-point func¬ 
tion, probably numerically on a lattice. This has been done for 
scalar field theory in 1+1 dimensions 2,3 , and more recently in 3+1 
dimensions within a large N expansion 4 . 

The 2PI approach is much trickier in a gauge theory than in a scalar 
theory, where most applications have been made. In a scalar theory it is 
generally assumed that the importance of a diagram, in the $ functional, is 
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determined by its loop order. But this is not completely obvious. And in a 
gauge theory, it is not true. The problem is that there are special kinematic 
ranges in a gauge theory, namely soft momenta and collinear momenta, 
where enhancements occur which can obviate naive loop counting. For 
instance, in QED, I claim that all the diagrams shown in Fig.l 





+ ... 


Figure 1. Some QED diagrams needed at leading order. 


are important at leading order. The proof that these are leading order is 
to open the photon line, to get a photon self-energy: 



Figure 2. Same diagrams, with photon line opened. 


and then to read our paper 6 , which does the detailed power counting. 



Figure 3. Left: scalar 2PI diagram which gives 2 <-> 2 scattering. Right: QED diagram 
which gives an interference term between scattering amplitudes, in more than one way. 


Cutting a 2PI diagram gives a matrix element and its conjugate. Al¬ 
ready knowing what scattering matrix elements we need for the Boltzmann 
equation, we can work backwards to figure out what 2PI diagrams must be 
required. For instance, a 3-loop diagram leads to interference effects which 
are important at leading order in QED, as shown in Fig. 3. 
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The trouble is that the inelastic, LPM suppressed splitting processes 
arise from the sum over an infinite set of diagrams; for instance, the inter¬ 
ference between two of the amplitudes for bremsstrahlung with 3 scatterings 

Cut 

\ 


A A A 

A A 

Figure 4. Two processes whose interference is important in determining the gluon 
bremsstrahlung rate, and the 2PI diagram they imply. Red lines are soft ( gT ) gauge 
propagators in the Landau cut. 




is shown in Fig. 4. In this and following figures, dotted lines are soft space- 
likc gauge boson propagators (in the Landau cut), while solid lines are 
always hard and on-shell; I use wavy rather than curly lines for gluons to 
make the pictures less cluttered. I have shown the loops responsible for the 
spectral weight of the soft lines for clarity; all lines in the 2PI formalism 
are always self-energy resummed. This is one of a family of diagrams which 
must be included. Any number of soft lines are permitted, provided that 
they do not cross; so for instance, in the calculation of shear viscosity it is 
necessary to include 2PI diagrams such as those shown in Fig. 5, just to get 
leading order results in g s - 



Figure 5. Examples of required diagrams. 


If one is willing to make extra approximations, the presence of this in¬ 
finite class of diagrams need not be fatal. There is a parametric difference 
between the soft exchange momenta and the hard momenta of the main 
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lines in the diagrams. The soft exchanges can be resummed. This leads to 
the integral equation for the splitting process already presented in Eq. (5); 
together with a gradient expansion, this should be the most efficient way 
to derive the Boltzmann equation presented in the last section. However, 
for the take-no-prisoners approach, an infinite set of diagrams is a show- 
stopper. Either some approximation must be made to resum them (and this 
probably demands some gradient expansion treatment or separation of time 
scales), the set of diagrams must be truncated, or the numerical approach 
becomes impossible. Truncating the set of diagrams raises new problems, 
because it is only the whole set, not any truncation, which is gauge invari¬ 
ant, and even then only when the distinction between gT and T momenta is 
made parametric. (A truncation also obviously means abandoning a leading 
order treatment-in fact, even a next to leading log treatment.) 

In conclusion, the 2PI formalism looks like a good way to derive the 
Boltzmann equation. It may provide a good framework for a consistent 
power counting effort, to prove that the claimed set of diagrams is complete. 
But as a direct computational tool, it appears ill suited in a gauge theory. 


3.2. Diagrammatic approach 

The diagrammatic approach to evaluating transport coefficients, in rela¬ 
tivistic field theory, was pioneered by Jeon 11 . Recently his (inefficient) 
treatment has been streamlined 12,13 and applied, at leading logarithmic or¬ 
der, to gauge theory 13 . Jeon showed that the diagrams which contribute at 
leading order are ladder diagrams, with a general form shown in Fig. 6. 



• Pinching 
, Propagators 

Rung 
(4-point 
_ function) 


Figure 6. General structure of ladder diagrams. Diagrams with any number of rungs 
(0,1,2,3,...) are needed. The crosses are operator or current) insertions. 


The self-energy needed on the pinching, nearly on-shell propagators is 
determined by differentiating the functional with respect to the propa- 
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gator once (opening one line), while the rung is obtained by differentiating 
twice (opening two lines). In scalar field theory, this leads to a relatively 
tidy set of diagrams. But in a gauge theory, the very complicated 2PI dia¬ 
grams discussed above lead to a set of diagrams, required at leading order, 
which is rather intimidating. For instance, in Fig. 7 I show the result of 



Figure 7. A bubble, which when opened gives a horrible ladder rung. 


opening one of the required 2PI bubble graphs on two of its lines, to get a 
required ladder rung. 

An illustrative example of one of the graphs in the general class which 
must be resummed, to determine the shear viscosity of QCD with fermions, 
is shown in Fig. 8. This diagram will hopefully give the impression that 




Hard off-shell 

Soft, spacelike, gauge boson, HTL resummed 
Hard on-shell, resummed with diagrams of form 



Figure 8. A typical diagram needed in the calculation of leading order shear viscosity 
in QCD. 
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computing transport at leading order in g s , in ultra-relativistic gauge the¬ 
ories, by a direct diagrammatic attack, is pretty hopeless. 

4. Conclusions 

Don’t get me wrong. I actually quite like the 2PI approach to dynamics. 
But it is not well suited to gauge theories, except as a method to do a 
convincing derivation of the Boltzmann equations, which should then be 
used to solve problems like transport coefficients. Direct numerical appli¬ 
cation on a lattice, such as has been done in scalar theories 2,3,4 , is probably 
hopeless (if complete leading order in g s results are desired-a Boltzmann- 
Vlasov approach, valid for transport at leading log order, may be feasible 
and might have some utility). 

As for the diagrammatic approach: I defy its supporters to give a con¬ 
vincing power counting argument, without reference either to 2PI or to 
kinetic theory, that the set of diagrams which I describe above, is both 
necessary and sufficient. Their resummation is expected to be technically 
challenging, and in the end it will only boil down to a derivation of the 
Boltzmann equations. I don’t see the benefit of doing this exercise. 
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